Electromagnetic Magic: The Relativistically Rotating Disk 
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A closed form analytic solution is found for the electromagnetic field of the charged uniformly 
rotating conducting disk for all values of the tip speed v up to c. For v = c it becomes the Magic 
field of the Kerr-Newman black hole with G set to zero. 

The field energy, field angular momentum and gyromagnetic ratio are calculated and compared 
with those of the electron. 

A new mathematical expression that sums products of 3 Legendre functions each of a different 
argument, is demonstrated. 

PACS numbers: 03.50.Dc, 02.03. Em, 02.30.Gp, 14.60.Cd 



I. INTRODUCTION 

Away from sources we write the electromagnetic field 
F = E + iB = -V* where * = <E> + iX is a complex po- 
tential satisfying V 2, F = 0. The simplest solution is the 
real one 'F = q/r but if we make q complex we introduce 
a magnetic monopole which is unacceptable. Moving the 
charge to a we have a potential qj \J {v — a.) 2 which cor- 
responds to a charge at a when a is real, but we now 
consider a to be complex ri + i&. Without loss of gener- 
ality we may move the origin to ri and orient the z axis 
along a. If we then set z = r cos 9 = r/i we find 



?/V(r 



q/(r 2 - 2iarn- a 2 ) 1/2 
1 \ 



(1) 



— Pnin) for r ^ a . 
r 



>F has a singular ring at fi = and r = a. To specify 
'F properly we need to say which square root is to be 
taken, but the expansion above for r ^ a is required to 
make everything regular at infinity and the simplest cut 
is across the circle defined by the ring singularity. With 
such a cut analytic continuation defines \& everywhere. 
The 'magic' field is given by 



F = E + iB = -V# = q 



ISl 



[(r - 2a) 2 ] 3 / 2 



(2) 



By construction V 2, F = except on the singular ring and 
the cut, so the sources of the field lie there. The expan- 
sion JJJ shows that the field has charge q, no magnetic 
monopole but a magnetic dipole qa. Indeed the electric 
field has only even multipoles and the magnetic field has 
only odd ones. On z = 0, a.r = so the denominator of 
P|l is real if r a and pure imaginary for r < a. 
Thus close to the plane of symmetry 



E = -qa/(a 2 - r 2 ) 3 / 2 



and B = —qr/(a 2 — r 2 ) 3 / 2 for z > and r < a. Below, 
where z < these fields are reversed. Notice that E is 
orthogonal to the disk so the disk is an equipotential and 



indeed its potential is zero (earthed) as may be seen by 
taking the real part of on z = with r 2 < a 2 . The 
magnetic field lies parallel to the radius below the disk 
and anti-parallel above. It does not cross the disk except 
at the singular ring. For r > a and z — 0, B points 
downwards everywhere. Thus every field line returns to 
the upper hemisphere through the singular ring. 

From the electric field one readily calculates the surface 
density of charge on the disk (summing both sides) 



qa 



27r(a 2 - r 2 ) 3 / 2 



(3) 



Notice that this is of the opposite sign to the total q 
which is rectified only at the singular ring. Indeed the 
total charge on the cut within r = R is Q{< R) = 

— q ( a ^ 2 a ^ fl r ~) for R < a which diverges with a negative 

sign as R approaches a but Q(< R) = +q for R > a so 
the negative infinity is now replaced by a finite positive 
result. Likewise the surface current in the cut is 



2tt (a 2 



2)3/2 



(4) 



which is precisely what we would get if the above charge 
density were rotating rigidly with angular velocity Q = 
c/a. The total current within R < a is Q(< R)Q/(2tvc). 
Again the magnetic effects of this current are over- 
whelmed by the current around the singular ring which 
is of opposite sign. Notice that the disk is not crossed 
by a magnetic field line saving at the singular edge ring 
itself. In this respect the disk acts like a superconductor 
displaying the Meissner effect. In the Black Hole context 
this phenomenon was noted by Bicak & Janis [I] (see also 
Bicak & Ledvinka [2]). 

I now list other properties of this 'Magic' electromag- 
netic field. (For proofs see Lynden-Bell [3]) 

1. Relativistic Invariants F 2 = E 2 - B 2 + 2iE.B = 
q 2 /(r — ia) 4 . 

2. E 2 = B 2 only on two spheres of radius v2a centred 
on r = ±a. They meet on the singular ring. 
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3. E.B = on the sphere r = a and also on the plane 
z = 0. 

4. The field energy density F.F*/8tt = [q 2 / (8ir)](r 2 + 
a 2 )/|(r — ia) 2 | 3 . This diverges like (r — a)~ 3 when 
z = near r = a. 

5. F* x F = 2igExB = 2iq 2 ax r/|(r-ia) 2 | 3 likewise 
diverges. 

6. The above is related to the Poynting vector which 
shows that the energy density flows around the axis 
with a velocity 1~2 x r where ft is constant on spheres 
r =const with fi(r) = pq^frr- 

7. Landau & Lifshitz [4] show that in the frame that 
moves with the velocity cV where V/(l + V 2 ) = 
E x B/(E 2 + B 2 ), the transformed fields E'&B' 
are parallel. Gair proved the theorem that these 
velocities are a uniform rotation of each spheroid 
confocal with the disk at the rate f2 = ac/ (f 2 + a 2 ) 
(for f see section 2). 

8. In the corotating frame of those spheroids E' and 
B' are both perpendicular to the spheroid on which 
they lie. 

9. The total field energy and the total angular momen- 
tum in the field both diverge due to their divergence 
at the singular ring. 

More remarkable properties still come from the sepa- 
rability of both wave equations and equations of motion 
of any relativistic charged particle in this field. 



answers for total field energy, total field angular momen- 
tum etc. We aim to find fields which give the Magic 
field as a limiting case but whose other members give fi- 
nite answers. Some spice is added to the investigation 
by Carter's [5] remark that all Kerr Metrics have the 
same gyromagnetic ratio, 2, as the Dirac electron, and 
Pfistcr & King [15] propose that this may have a deeper 
significance for general relativity. Indeed it holds for all 
the conformastationary metrics, the statement by us [16] 
that the gyromagnetic ratio was one, was based on erro- 
neously missing out the factor 2 in its definition. Our disk 
[16] remains the only known conformastationary interior 
solution to Einstein's equation. 

Recently [17] we discussed the remarkable behaviour of 
charges on a relativistically rotating conducting sphere. 
The changes in the field correspond to addition of more 
and more of the Magic field as the rotation increases. 
However even when v = c the distribution on the sphere 
never becomes the Magic field: that discussion leads to 
the suggestion that the Magic field may be the field of 
the rapidly rotating charged conducting plate with a tip 
speed of c. The facts that the component of E + V x B 
in the plate is zero and that the current is due to the 
convected charge strongly suggest that this is the limit, 
as Vc the tip speed tends to c, of the charged rotating 
conducting plate. This paper is devoted to a discussion 
of the whole sequence of fields of such a plate when it 
rotates at any tipspeed up to c. The field lines of the 
Magic field are shown in Figure 1, which was generated 
by J. Gair. 



1. The Hamilton- Jacobi equation separates 

2. The Klein Gordon equation separates 

3. The Dirac equation separates 

4. The Schrodingcr equation separates if only the real 
part of the field is included (but not when the mag- 
netic part is added although the reverse is true for 
the Klein Gordon equation). 

All the above really stem from remarkable investiga- 
tions in the separability of wave equations around black 
holes by Carter [5], [6], Teukolsky [7], Chandrasekhar [8] 
and Page [9]. In the simpler flat space case see Lynden- 
Bell Paper I [10]. 

Finally the relationship with black holes and indeed the 
original discovery of the field in General Relativity is due 
to Newman et al. [11] who generalised the Kerr metric of 
a spinning black hole to include charge. If in his solution 
one puts Newton's G = one obtains an electromagnetic 
field in flat space which is the above 'Magic' field [12], 
[13], [14]. 

The present investigation is aimed at providing a set of 
electromagnetic fields in fiat space that keep some of the 
remarkable properties of the Magic field but give finite 



II. MAXWELL'S EQUATIONS FOR THE 
ROTATING DISK 



We treat the problem using oblate spheroidal coordi- 
nates confocal with the ring that forms the edge of the 
disk. Actually we shall solve the problem as the limit 
of a highly flattened spheroid as its minor axis b tends 
to zero; this allows us to make direct comparisons with 
our treatment of the rotating charged conducting sphere. 
The equation 



= 1 



(5) 



defines a set of confocal oblate spheroids each with a 
definite positive value of r. When r = it defines a disk 



0. 



x 2 + y 2 



^ a 2 . When r is negative it defines a set 
of hyperboloids of revolution (of one sheet) each of which 
is orthogonal to the oblate spheroids. For r positive we 
write it as r 2 and for r negative we write it as —a 2 cos 2 ■& 
or more often as —a 2 fi 2 . The spheroidal coordinates r, d 
are supplemented by the usual azimuthal angle (f> to give 
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FIG. 1: (a) Electric lines of force for the Kerr-Newman potential. The density of lines reflects the strength of the field, 
(b) Magnetic lines of force for the Kerr-Newman potential. The density of lines reflects the strength of the field. 



the metric where 9 is the spherical polar coordinate which should be 

^2 _|_ fl 2 ~2 i~2 _|_ fl 2 ~2\ distinguished from the spheroidal coordinate i? although 

ds 2 — — — — a ^ dr 2 + — a J^ dp 2 they are equal at oo. We consider the surface of the 

r + a 1 — M spheroid r = b to be the rotating conductor. 

+ (1 — A 2 )(^ 2 + a 2 )d(f) 2 (6) On f = 0, the disk, dX/df changes sign which corre- 

= dr 2 + r 2 d6 2 + r 2 sin 2 6d<j) 2 sponds to the continuity of the normal magnetic flux, and 

= dR 2 + R 2 d(j) 2 + dz 2 
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$ is continuous which corresponds to the surface compo- 
nent of the electric field being continuous. 

The other requirements are that the total charge must 
be q, so $ — > q/r at oo and E + i(l~2 x r) x B must have no 
component along the disk's surface so that the Lorentz 
force component vanishes. Of course the discontinuity of 
the surface component of — VX corresponds to the sur- 
face current AttJ^ and the discontinuity in the normal 
component of — V$ gives the surface density of charge 
Ana. In the steady state with an imperfect conductor 
the currents are just those caused by the rotation of the 
charge. 



J = — (SI x r)cr 

c 



(7) 



That completes the physical specification of the prob- 
lem. Now we turn these physical requirements into math- 
ematics. 

On the disk 

Jm(-d^/dR) = B R = 2ttJ 
TZe(-d^/dz) = 2ira 
and from the Lorentz force component 

Tle(-d^/dR) + (SlR/c)Xm{-d^ /dz) = 
while 10 gives 

lm(-d^/dR) - {SlR/c)TLe{-d^/dz) = 
Hence on the disk R < a, z = 

d^/dR=i{SlR/c)d^/dz (8) 

Now in spheroidal coordinates the general solution of 
Laplace's equation which tends to q/r at oo is 



la ' 



(9) 



Pi{fj) are the Legendre Polynomicals and the Qi (■?-) are 
the Legendre functions here expressed as functions of an 
imaginary argument. The Q21 (£•) are imaginary while 
the Q21+1 (-?-) are real. The Qi(z) obey the recurrence 
relations (/ + l)Qi + \(z) = (21 + l)zQi(z) — IQi-i(z) and 
(z 2 -l)Q' l (z) = lzQ l (z)-lQ l _ 1 (z). 

Hereafter whenever we refer to Qi or Q[ without ex- 
plicit mention of its argument we shall mean Qi (0) and 

Q[(0) i.e. the limit of Qi (i) as b -> etc. We then 
have 

Ql-i = -[(l + l)/l]Qi+i and Q\ = IQi-x = -(l + l)Qi+i 

(10) 



In 



h(z) = -\n 



- 1 



in b 

as z — > > 

2 ia 

■ 1- f 

la/r as z = > 00 



— 1 as z 



The condition that "J - — > q/r at 00 becomes \&o = 1. 
Rewriting our boundary condition (JSJ in spheroidal co- 
ordinates 



v/1 -ji 2 _ ifia v/l-/2 2 5* 



a|/i| 9/2 c 
Thus from setting Qa/c = u> 



Or 



on f = 



Now as well as their usual orthogonality relation the 
Legendre polynomicals Pi(fi) have the property 



2 for I ^ m, (I + to) even 
for I < m or (7 + m) odd 



Multiplying (fTT|) by |P2m(A) and integrating —1 to +1 



^(4Z+3)?72i+i <32i+i = -wr? 2m Q 2m = -wife ro 2mQ 2 m-i 
m 

(12) 

where 77 m = f m / (2m + 1). Subtracting this with (m + 1) 
written for m 

(4m + 3)772m+lQ2m+l = W [7?2m+2 (2m + 2)Q 2 m+l 

+ V2m(2m + 1)Q 2m+lJ 

where we used (|10|l on the last term. Thus 

(4m + 3)?72m+i = w [(2m + 2)77 2m+ 2 + (2m + 1)772™] 

(13) 

similarly multiplying l|ll|) by iP2m+i(M) an( i integrating 
etc 

(4m + 5)?72m+2 = w [(2m + 3)772,^+3 + (2m + 2)ryam+i] 

(14) 

(|13|l & (|14ll can be recombined into a single formula true 
for all m 



(2m + l)r) m = lu [(m + l)?/ m +i + mr/ m _i] 



(15) 



When the Magic potential is rewritten in terms of the 
spheroidal coordinates f and /i it takes a remarkably sim- 
ple form and it is interesting to resolve it into Legendre 
polynomials of jj, 



^ = q/y/ (r — ia) 2 = q/(f — iap) 



We multiply both sides by P n (P>) and integrate —1 to +1. 
By Abramowitz & Stegan 8.8.3 \ ^j^y^A 1 = Qn(z) 
so we find Mi (j 1 ) = (21 + l)Qi (-?-) and in terms of the 
expansion @ = 21 + 1 so r\i = 1 for the Magic field. 
We see at once that this satisfies <|15[) with ui set equal 
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to 1 so that the Magic held is as expected the limit of 
the uniformly rotating disk when v — ► c. To solve our 
problem for a general rotation rate we use the generating 
function method developed in the appendix of Paper III 
[17]. 

We multiply equation 1)15(1 by U m and sum from m = 1 

oc 

to oo. Defining r](U) = ^2r) m U m and remembering that 
o 

7/0 = $o = 1 we find that 

2Urj' + rj - 1 = u)(r)' - 7)i + U 2 rj + rj) 



We rewrite this in terms of the roots u ^ 1 and 1 ju of 
the quadratic t 2 - 1 = so ^ = ^(u^ 1 +u). [Notice 
that the quadratic and u differ from those of paper III in 
which we discussed the sphere rather than the disk] 



(U -u)(U- u- 1 )t)' + [U- liu- 1 + u)]r) = T)i - - (16) 

2 u 

The general solution for the generating function rj(U) is 



i^Km--) / r'du + c] 



(17) 



where /, the integrating factor, is given for U ^ u by 
I=[(U- u){U - u- 1 ) 1 / 2 = [1 - -U+ U 2 ] 1 ' 2 . (18) 



In general the expression (|17fl has a factor be- 
having like (it — t/) -1 / 2 which diverges at U = 
u. Its expansion as a power series is of the 

form it- 1 / 2 ( 1 + ku^U... + i/ii%iSz( u - 1 [/)™ + 



2" ' n\ 

which gives rise to terms in it~( n+1 / 2 ) which diverge as 
n becomes large because u < 1. To suppress this diver- 
gence we must choose the constant 



C = 



i)i - - 

to 



r l du 



Then the numerator vanishes at U = u and 

1 



V : 



But 770 = 77(0) = 1 so 771 



r 1 / i- x du 
-[f u i-i d u}-\ 



Now I 



1(0) = 1, so we set ^(u^ 1 + u) - U 
then / = ^(u^ 1 — u)sh w and w — when U = u. Also 



u)] and 
i(it _1 — u)ch w; 



r x du 



w 



Hence 



and 



T){U) = - 



l(u~ 1 + u) 



w 

shw 



(19) 



(20) 



where 7 2 = 1 — lo 2 = \ — v 2 /c 2 ; 7 
77(0) = 770 = 1 so we have 

r)i — 1/u) — XjchT 1 ^ — > 1 as ui - 



But 



(21) 



Thus at u> = 1 our recurrence relation H15[l yields the 
Magic field with r\ n = 1 all n. We now know both 770 = 
1 and 771(0;) as a function of u>; thus we are now in a 
position to use our recurrence relation to generate all the 
r/ m for any chosen value of to. Alternatively we could 
expand our expression for rj(U) as a power series in U 
and pick out the coefficient of U m . However Prof J.F. 
Harper pointed out to me that the recurrence relation 
((15JI is that obeyed by the Legendre Polynomials P n (^). 
Notice that ~ is greater than 1. Since the P n and the Q n 
obey the same recurrence relation the general solution is 
a linear combination of them. Fitting that combination 
to 770 = 1 and the value of 771 just derived we find, noting 
thatc/i- 1 7 =Qo(^), 



'7 = <?„ ( — J /Qo ( — 



Thus the solution for the complex potential is 



(22) 



Qr, 



LO I \LO 



Q, 



~ Pn(A) 

ia ' 



(23) 



We note that Q n (^) /Qo (~) — > 1 as lo — > 1 so the 
Magic field is still there in that limit. Now by [18] 8.8.2 
Qn(z) —> z~ n ~ x f Q (1 + cht)~ n ~ 1 dt as z — > 00. Hence 



Q, 



i) /Ob (i 



\2n+2 



so for small w 
* = ■ 

Now 

Qo (- 

V ia 



Qi(- 

70 



/I (1+^) 2 

r Lo n 2 n (n\) 2 /(2n+ 1)! 
+0(w"+ 2 ) 



(2n)! 



(24) 



In 



i tan 



-1 (t 

V f 



Q2 



— tan — — 1 ; 

a f 

7.2 l„2 



-3f 2 /a 2 



i tan 



-(f) 



z3r 
2^ 



so 




ft (A) 



o(^ 3 )} 



(25) 
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The first term is the well known potential of a static 
charged disk. 

We may now read off expressions for the charge q, 
the (magnetic) dipole moment fi m and the (electric) 
quadrupole moment Q e 



\qaw + 0{u 2 ) 
qa ; lo — ► 1 




where we have defined X = In «/ = ch 1 -f. 

In finding the quadrupole moment for u) small we note 
that part of it comes from the oblateness of the f =const 
surfaces. Taking the sign so an oblate distribution has 
negative Q 



Qe = -^a 2 



I +-2Q-2 ( i ) /A 



\qa 2 (X 

,2 



5 W J 



-gar 



III. CLOSED FORM EXACT POTENTIAL 

On the axis p, — 1, so all the P n (p) — 1 and we can then 
perform the summation (12.'i[l using a result in Whittaker 
& Watson [19] 



l)Qr. 



2(^-i) 



r 
ia 

X + u f — ia 



X — u) f + ia 



On the axis f — z so we obtain the complex potential 
setting b = — 



9{z) 



ib 



1 , / z — ia 
1 + — In 

2A V z + ia 



(26) 



Now there is a remarkable way [20] of obtaining the so- 
lution of V 2 ^ = everywhere if it is given as an analytic 
function ^f(z) on the axis. The potential everywhere is 
given by (see Appendix A) 

1 f 2 * 

W(R,z) = — / V(z + iR cos a)da . 
2tt Jo 

In our application it is simplest to use tan ^ = t as the 
variable of integration. The solution is then 

, , 1 f°° ( l-t 2 \ dt 

applied to the ^>(z) given by (|2*S|) we find 

.(*.)-!/" 



1 + 2a'° 



(z — ib+ iR) + (z- ib - iR)t 2 
(z — ia + iR) + (z — ia — iR)t 2 



(z + ia + iR) + (z + ia — iR)t 2 



dt .(27) 



Now 



t 2 



fJ 



and 



ln(/ 2 +t 2 ) 



t 2 



dt 



9 



9) 



(28) 



(29) 



the latter result may be checked by differentiating with 
respect to f 2 performing the integral that results by par- 
tial fractions and finally reintegrating with respect to f 2 . 
The constant of integration is readily seen to be correct 
by taking the limit f g; for z > and R > no 
singularities occur in the t integrations. We are most 
interested in the region R ^ a ^ b and z small. We write 



and 



fl 
I'- 
ll 2 = 



z + iR — 


ib 


b 


-R + 


iz 


z — iR — 


ib 


b 


±R + 


■ 

iz 


z + iR- 


ia 


a 


-R + 


iz 


z — iR — 


ia 


a 


+ R + 


iz 


z + iR + 


ia 


a 


+ R- 


iz 


z-iR + 


ia 


a 


-R- 


iz 


a + R + iz 






- iR- 


- ia 


a — R — iz 




z 


- iR- 


- ia 



Writing z = rcosO, R = rs'm9 we see that for large 
r, g, f + and /_ all become e l6 while h is e ITr / 2 . On 
the axis R — 0, g, /+ and /_ are all unity while h 

T / 2 



is e 



which becomes 1 as 



z +ia — — > 0. Performing 

the integrations and making sure that we take the right 
branch at infinity we find 



4- = 



<1 



X 



ln(-ih) + In 



(z — ib — iR)g 
so we have the solution to our problem 
<1 



f± + .9 

f-+g 



ib) 2 



l 



l 



ITT 

2A 



Defining S(a, R, z) — \J a + R + iz 

S{a, -R, z)S{b, R, z)+S{a, R, z)S{b, -R, z 



(30) 



Z = ln 



S(a, R, —z)S(b, R, z) 



-S(a, —R, —z)S(b, —R, z) 

(31) 



Although this expression looks a bit formidable its deriva- 
tives at z = where we shall need to evaluate them are 
quite simple 



0Z_ 

dz 

dZ_ 
OR 



ib 



R 2 ){b 2 
R 



R 2 ) 



o y/(a 2 -R 2 )(b 2 -R 2 ) 



(32) 



(33) 
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We have written the forms appropriate for R a. 

In expressions (I3U|1 and l|31[) we have demonstrated how 
the formidable sum (|23|l is evaluated explicitly but we 
are not aware that this identity has been demonstrated 
by experts in mathematical functions. 

It is now simple to check that the boundary condition 
which may be rewritten 



0* 
dR 



iR a* 

b dz 



for R ^ 



is indeed satisfied by our solution H3()|l . One consequence 
of this is that it is easy to integrate the surface density 
of charge 

a= ne [- — 

ZTT \ OZ j q 

evidently 

rR r R b a* 

2naRdR = Re / — ^dR 
h, Jo idR 

= -6Xm[*(i?,0) - *(0,0)] 



Q« R) 



<? 1 



1 



Z(R,0) 
A 



y/b 2 - R 2 

We note that Z(a,0) = —A because b = a/u> 

expected all charge lies at R ^ a. As ui — > 1 

In \ 1 - for all R < a but A — > oo 



so as 



lim 



f Z(R,0) 



A 



for i? < a 
-1 for R = a 



which explains the somewhat strange behaviour of Q(< 
R) for the Magic field. 

At the pole the surface density of charge becomes zero 
when A = -j which occurs at ui = 0.834 which gives a tip 
speed of 83.4% of c. This should be compared with 93% 
for a sphere. The difference is in the expected direction 
since the disk has an uneven charge when lo — and for 
it the Lorentz magnetic force lies radially across the disk. 



IV. FIELD ENERGY, ANGULAR MOMENTUM 
& GYROMAGNETIC RATIO 



by 



The total energy in the electromagnetic field is given 



e = (Sn)- 1 J (E 2 + B 2 )dV 

= -(Sir)- 1 [ ^*d^/dz4nRdR 
Jo 

= i(b/2) [ V*dV/dRdR 
Jo 

= -6/2 Im / ®*Tr=:dR 
Jo dR 



Now 



9* 
dR 



qR 



3/2 



{b 2 - R 2 ) 



7T ./ Z 

2A + U 1+ A 



iqR 



AVa 2 - R 2 (b 2 - R 2 ) 



and 



4r* 



q 



Vb 2 - R 2 



2A 



i 1 + 



so apart from the last term in d^f /dR, we see that ty* and 
d^ /dR have opposite phases. Thus only the last term in 
d^ /dR contributes to the imaginary part of ^>*d^ /dR. 
Hence 



irb 
8A2 



2RdR 



o (a 2 - R 2 y/ 2 (b 2 - R 2 ) 3 / 2 



The integral is 2(a/6)(6 2 - a 2 ) -1 so 



7r aq 

4A2 (62 _ a 2 



= 7 r( (? 2 /a)- 



(1-^) 



1 1— u; 



(34) 



which gives the correct limit of j(q 2 /a) when u> — > 0. 

The total angular momentum in the electromagnetic 
field is 

L = — / rx(Ex B)dV = — !— / r x (V** x V*)dK 
Aire J 8iric J 

We note that the integral has no real part. 

The Angular Momentum can be rewritten in the form 

L = — / rx(Ex B)d 3 x 
Aire J 

= Y~Zm try. [W* x (VAR x V0)]d 3 x 



-Im / r x V^Vr.VfARW 3 

7TC 7 

Im J W*.V{AR)d 3 x 



8ttc 

z 
8ttc 



but since V 2 ** = 0, div(ARV**) = V**.V(AR) so we 
may turn L into a surface integral 



L = — Im / ARdVVdzAirRdR 
8ttc 7 

where the factor 2 accounts for both sides of the disk. 
Now on the disk Rd^*/dz = ibd^*/dR by © so on 
z = 0, 
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—Tie 

2c 



zb. 
~2c 
zb 

Yc 

zb 
2c 



Tie 



Tie 



ARd^*/dRdR 



-**(a, o)A(a, o)a 



^*d/dR{AR)dR 



^*Aa 



<9* 
*£*R—dR 

az 



ne[-aV*A] a -lm[b J **<9*/<9i?di?] J 

We already evaluated — Im b f° ^* / dRdR in con- 
nection with the energy it gave 27r(g 2 /a)- 



soL = nz(q 2 /c) 



(x^^-^^a-Toeval- 

uate the last term we need to find A from the formula 
12|) of the Appendix putting tana/2 = t this becomes 



A(R,z) = ^ 

TT 



(1 + t 2 )[(z - it + iR) + ( 



(z — ia + iR) + [z - 



z 
ia - 



ib - iR)t 2 } 
iR)t 2 



(z + ia + iR) + (z + ia — iR)t 2 
Evidently we need the following integrals 
1 f°° 1 - t 2 ln(/ 2 + t 2 ) 



dt 



I = 



1 



t 2 
■In 



g 2 + t 2 
/ + 1 



-dt 



(g-1) 

(ff 2 -l) i "V/ + .9; .9(5 + 1) 



Mf + g) 



This expression is derived by first differentiating / with 
respect to f 2 then integrating over t and finally integrat- 
ing with respect to f 2 . The final 'constant' of integration 
which might depend on g, is found to be zero by looking 
at the full limiting values when / 3> <?. 



K = 



1 



1-t 2 



1 



1 + t 2 g 2 + t 2 



dt = - 



g-1 
0(0 + 1) 



Using these forms to evaluate A we find 



A(R,z) 



iq 



(z — ib — iR) \ g(g 
In h + In ' ' 
2 



1)A 



ITT 

~2 



f-+9 

/+ + 1 



(9 2 ~ 1)A 



In 



f- 



1 



where / + , /_, g and /i have the same definitions as be- 
fore, and Inh + ln[(/+ + «?)/(/_ + <?)] = Z. 

Although the above expression gives the complex A ev- 
erywhere, we only need A at (a, o) to evaluate the angular 

b— a / 1— uj ^ 

b+a ~~ V 1+w ' " 



moment um . There Z = 
oo like 



-A, g 



and /_ 
final In term in A becomes ln(l/g) 



-=%. With these substitutions the 

a — it 



A. 



Evaluating A(a, o) we find 



A = 



VT 



Z7T 

-1 — +1 



Now from ((301) *(a,o) 



and 



^e(A«-* 



Trquj 



Hence 



ab 



nq LU 

2^yr=^A 



1 



2c'~ v ~~~ ' 4c VI - w 2 A 

I remind those who are worried that this may appear to 

have the wrong behaviour with to that A = In <J is 

odd in lu and ~ lu for lu small. Putting this result into 
our expression for L 



L = 7TZ l(oj) 

C 



(35) 



where 



uj 2 In ( tta 



(1-w 2 ) ln( 



1— w 



This expression tends to L = 7r^§^ for small lu. The 
dimensionless quantity that is the inverse fine structure 
constant in the case of the electron is 



2Lc/g 2 = 2ttI(lu) 



(36) 



This assumes that all the angular momentum is electro- 
magnetic. An lu of 0.999743 gives the right value for the 
fine structure constant; this corresponds to an equatorial 
Lorentz factor 7 of 44.1. 

The gyromagnetic ratio is [A m /L and the g m factor is 



2mc Mj? 



9 m = 2 



l-Vl-w 3 A/u 



with A 



iln 



l + UJ 

1-u 



(37) 



As 



where G(u>) = 

lu — > 1, A — > 00 but vl — cj 2 A — > 0, so G — > 1 and 
the gyromagnetic ratio would tend to 2 if all the energy 
were electromagnetic. However there is an unbalanced 
electromagnetic stress on the disk's edge and most ways 
of making a balancing stress also add to the energy. Some 
add to the angular momentum too. Perhaps the simplest 
way to a consistent model is to add a string loop around 
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the edge of the disk with energy per unit length equal 
to its tension T. Such a string loop carries no angular 
momentum but has energy 2iraT. The total energy of the 
whole configuration is then e + 2iraT. Minimising this 
over 'a' while keeping the angular momentum L fixed 
via (|36|) means keeping w fixed. Hence from Ij34|l e is 
proportional to a -1 and equilibrium at the minimum is 
achieved when 2nTa = e. Thus the total energy is 2e 
and mass of the whole configuration at equilibrium will 
be 2e/c 2 ; for such a model the gyromagnetic ratio or 
rather the g m factor is 

g m = 4G(w) 

which becomes not 2 but 3.46 for the w that gives the 
correct fine structure constant. Probably it is more sense 
to model the disk with a stress-energy tensor based on 
the Kerr disk (Lasenby et al. 2004 [21]) which contains 
both angular momentum and energy. 

It is sensible to remark that infinities such as that 



found here as v — > c are often removed when relativis- 
tic problems are treated quantum mechanically. 

In [3] we gave singular solutions that rotated uniformly 
at any rate and were superpositions of forward and back- 
wardly rotating Magic fields. The equal superposition 
was static but did not give the well known solution of the 
charged static disk. This gave us doubts as to whether 
the Magic field would be the v — > c limit of the rotating 
disk problem but we have now shown it to be so. The 
other singular solutions are discussed in Appendix B. 
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APPENDIX A: SCALAR AND VECTOR 
POTENTIALS FOR POLOIDAL FIELDS 

We consider fields which lie in the meridian planes 
through an axis of symmetry and which are harmonic 
above the plane z = 0. For simplicity of exposition we 
take their potentials to be O(-) at oo but the method is 
readily extendable beyond those cases. Since the fields 
are harmonic we may write 



E- 



iB = -W = Curl A 



It is normally not hard to derive an expression for 'J/ on 
the axis of symmetry and we shall suppose this has been 
done so that ^(z) is known, and *$>(z) — » q/z as z — > oo. 

We show that the complex scalar and vector potentials 
^(R, z) and A(j> where <j> is the unit toroidal vector are 
given by 



1 f 2 " 

®(R,z) = — / *(z + iRcosa)da 
2tt Jo 



A(R,z) = — / cos a$(z + iRcos a)da 

provided ^(z) is analytic in z < 0. 
Proof: 



(Al) 



(A2) 



K ' 2m J 



<K 
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where the path in the complex £ plane is chosen to en- 
circle £ = z. 
Hence 

1 f *(0 

*(z -Hi? cos a) = — <p d( 

2m J Q — z — iR cos a 

and 

[yr(z + iRcosa)da = — [ I- d(da 

J 2iri J J Q — z — iRcosa 

We choose the integration circle to surround the whole 
circle z + iRcosa for all a. We write tana/2 = t and 
reverse the order or integration 




(C - z - iR) + (C - z + iR)t 2 " ^(C-z) 2 + i? 2 



7T 

= v/(r - Cf ) 2 

Thus 

1 /* 27r 1 /" *(0 

— / V(z - iRcosa)da = - — <p ■ ' - 
2tt J V ; 27tz J ^/(r-Cf) 2 

But V 2 , 1 is zero everywhere except at r = tr 
which is not encountered on the integration path. Hence 
5sr fo^ ^'( z + *-R C0S CK)(£r is harmonic. But it also reduces 
to ^>(z) on axis and is O (i) at oo so it gives the complex 
potential everywhere z ^ 0. 

To obtain the vector potential we note that 

CurlA = V(AR) x V<p = -V* (A3) 

and 

Curl CurlA = -[V 2 (AR) - 2V lniZV(AR)]V</> = 

for z > because there are no sources above z = 0. 

We now write A = dQ/dR then the expression in 
square brackets above reduces to 

i? d/dR(W 2 C) = 

so 

v 2 C = /(*) 

Now we may add an arbitrary function of z to £ with- 
out changing A so we add — JJ Z f(z)dz and deduce that 
the new £ will obey 

V 2 C - 



from (|A3|) and the above 

d 2 ( _ l_d_ [ R dC\ __ 9* 
~ dz 2 ~ R dR \ dRj ~~ dz 

Hence ((z) = J z ^{z')dz'. 

The linear function of z that comes from the integra- 
tion can again be absorbed into £ without changing A or 
the fact that V 2 C = 0. 

We now use the theorem proved above for taking ^ off 
axis, but apply it to Then 

C(R,z) = — J J ^(z')dz'da 

differentiating with respect to R 
i f 27T 

A = — / cos a^> (z + iR cos a)dxa . 
2tt Jo 

We notice that both £ and A obey their appropriate 
equations and that A is O (^) at oo as required. 



APPENDIX B: SINGULAR SOLUTIONS 

In [3] we remarked that if * = q[R 2 + (z - ia) 2 ]~ 1/2 
is the complex potential of the Magic field itself then 
a^> + (1 — a) 1 ^* is a superposition of a forward rotat- 
ing and a backward rotating Magic field. The electrical 
potential of the disk is still zero and no magnetic field 
crosses it in R < a. There is the same charge density on 
the disk as in W and it rotates uniformly but less fast. 
When a — i the charge does not rotate but the surface 
density is still singular. This is not the solution for the 
non-rotating charged disk but to what physical problem 
does it correspond? I believe it must correspond to the 
field of an earthed plate surrounded by a charged wire 
with a small insulating gap, between the wire and the 
plate. As the gap becomes smaller, a larger and larger 
charge must be placed on the wire to ensure that the net 
charge of wire and plate remains at a fixed value q. I 
think this non-rotating magic field corresponds to such a 
singular limit. Perhaps the rotating versions with | 
correspond to the singular solutions of the recurrence re- 
lations that we rejected. These solutions obey the con- 
ditions required of a superconducting disk, the charges 
while rotating uniformly do not correspond to the phys- 
ical problem that we set out to solve. Nevertheless the 
Magic field itself is the b-»c limit of our problem. 



